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Abstract 

We present simulations of evacuation processes using a recently introduced cellular 
automaton model for pedestrian dynamics. This model applies a bionics approach to 
describe the interaction between the pedestrians using ideas from chemotaxis. Here 
we study a rather simple situation, namely the evacuation from a large room with 
one or two doors. It is shown that the variation of the model parameters allows 
to describe different types of behaviour, from regular to panic. We find a non- 
monotonic dependence of the evacuation times on the coupling constants. These 
times depend on the strength of the herding behaviour, with minimal evacuation 
times for some intermediate values of the couplings, i.e. a proper combination of 
herding and use of knowledge about the shortest way to the exit. 
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1 Introduction 



Methods from physics have been successfully used for the investigation of ve- 
hicular traffic for a long time [1,2]. On the other hand, pedestrian dynamics has 
not been studied as extensively [3]. Due to its generically two-dimensional na- 
ture, pedestrian motion is more difficult to describe in terms of simple models. 
However, many interesting collective effects and self-organisation phenomena 
have been observed (see [2,4] for an overview and a comprehensive list of 
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references), e.g. jamming and clogging, lane formation and oscillations at bot- 
tlenecks in counterflow or collective patterns of motion at intersections. These 
phenomena will be discussed in sec. 2.4. 

The model takes its inspiration^] from the process of chemotaxis (see [5] for a 
review) . Some insects create a chemical trace to guide other individuals to food 
places. This is also the central idea of the active-walker models used for the 
simulation of trail formation. In the approach of [6] the pedestrians also create 
a trace. In contrast to trail formation and chemotaxis, however, this trace is 
only virtual although one could assume that it corresponds to some abstract 
representation of the path in the mind of the pedestrians. Its main purpose is 
to transform effects of long-ranged interactions (e.g. following people walking 
some distance ahead) into a local interaction (with the "trace"). This allows 
for a much more efficient simulation on a computer. 

The basic idea of our approach might be used for studying a variety of prob- 
lems, especially from biology [6-8]. Here we want to apply this model to a sim- 
ple evacuation process with people trying to escape from a large room. Such a 
situation can lead to a panic where individuals apparently act irrationally. A 
nice discussion of empirical results can be found in [4]. Our motivation here is 
rather the determination and classification of the different types of behaviour 
exhibited by the model than a realistic application. 

The phenomena observed during panics can be quite different from those found 
in "normal" situations. Nevertheless it is desirable to have a model which is 
able to describe the whole spectrum of possible pedestrian behaviour in a 
unified way. So far mainly the social-force model [9] has been used which 
allows to reproduce the observed behaviour [2,4,10] quite accurately. In this 
continuum model the pedestrians interact by a repulsive (social) force which 
decays exponentially with the distance between them. This means that in 
each step of a simulation of N individuals 0(N 2 ) interaction terms have to be 
evaluated. Furthermore, in complex geometries it occurs quite frequently that 
two pedestrians are rather close to each other but do not interact since they 
are separated by a wall. Therefore in principle one has to check for all pairs 
of individuals whether an interaction is possible or not. For large crowds this 
becomes very time consuming. In contrast, in the model used here pedestrians 
only interact with the floor field in their immediate neighbourhood. Therefore 
one has only O(N) interaction terms. A further advantage is the discreteness 
of the model which allows for a very efficient implementation for large-scale 
computer simulations. 

We start with a short summary of the models basic concepts. 



3 Such "learning from nature" is the central idea of a field called Bionics. 
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2 Basic principles of the model 



In the model the space is discretised into small cells which can either be empty 
or occupied by exactly one pedestrian. Each of these pedestrians can move to 
one of its unoccupied neighbour cells at each discrete time step t — > t + 1 
according to certain transition probabilities (see fig. 1). The probabilities are 
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Fig. 1. Allowed motions and the corresponding transition probabilities. 

given by the interaction with two floor fields [6] . These two fields 5* and D 
determine the transition probabilities in such a way that a particle movement 
is more likely in direction of higher fields. They will be defined in the following 
subsections where also the basic update rules of the model are summarised. 

2. 1 The static floor field S 

The static floor field S does not evolve with time and is not changed by the 
presence of the pedestrians. Such a field can be used to specify regions of space 
which are more attractive, e.g. an emergency exit or shop windows. In case 
of the evacuation processes considered here, the static floor field describes the 
shortest distance to a an exit door, lying at the middle of the top wall of the 
room. Fig. 2 shows graphical representations of 5* for different geometries. S is 




Fig. 2. Static floor field S: a) for a lattice of size X x Y = 33 x 33 with one exit of 
width five cells; b) lattice of size X x Y = 30 x 60 and four exits. 

calculated due to a certain distance metric for each lattice site so that the field 
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values are increased in the direction to the door. The field values are highest 
for the door cells. The explicit construction of S can be found in Appendix A. 



2.2 The dynamic floor field D 

The dynamic floor field D is a virtual trace left by the pedestrians and has its 
own dynamics through diffusion and decay. It is used to model an attractive 
interaction between the particles. Fig. 3 shows 3-dimensional plots of D for 
three stages during evacuation process. At t — for all sites of the lattice 




Fig. 3. Dynamic floor field D: a) beginning evacuation process; b) middle stages; 
(c) final stage of evacuation. 

the dynamic field is zero, i.e. Dij = 0. Whenever a particle jumps from site 
to one of the neighbouring cells, D of the starting place is increased by 
one: — > + 1. 

Therefore D has only non negative integer values and can be compared to a 
bosonic fiele^\ i.e. the bosons dropped by the pedestrians during their move- 
ment create the virtual trace. Thus the field value corresponds to 
bosons. The dynamic floor field is time dependent, it has diffusion and decay 
controlled by two parameters a G [0, 1] and 5 G [0, 1], which means broadening 
and dilution of the trace. In each time step of the simulation each single bo- 
son of the whole dynamic field D decays with probability 5 and diffuses with 
probability a to one of its neighbouring cells. Finally this yields D = D(t, 5, a). 

Note that the dynamic floor field D is only altered by moving particles and 
therefore it corresponds to a virtual velocity density field, rather then a particle 
density field. In sec. 3.2 we discuss an alternative definition of the dynamic floor 
field D as a pure particle density field and its consequences. In the following 
we give a short summary of the specific update rules used in the simulations. 



4 A simplified one-dimensional model with a fermionic field has been used recently 
for the description of chemotaxis [8] . 
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2.3 Update rules 



The update rules of the full model including the interaction with the floor 
fields have the following structure [6,11,12]: 



(1) The dynamic floor field D is modified according to its diffusion and decay 
rules (for details, see sec. 2.2). 

(2) For each pedestrian, the transition probabilities p^ for a move to an 
unoccupied neighbour cell (see fig. 1) are determined by the local 
dynamics and the two floor fields. The values of the fields D and S are 
weighted with two sensitivity parameters k$ € [0, oo[ and ko £ [0,oo[. 
This yields 

p i:j = Nexp (k D Dij) exp (k s Sij)(l - n,j)0.j , (1) 

with 



occupation number: 
obstacle number: 



riij = 0, 1, 



for forbidden cells, e.g. walls 

1 else 



normalisation: N 



ex P (koDij) exp (k s Sij)(l - »; ; )C;, 



(3) Each pedestrian chooses a target cell based on the transition probabilities 
Pij determined in the previous step. 

(4) The conflicts arising by any two or more pedestrians attempting to move 
to the same target cell are resolved by a probabilistic method^. The 
pedestrians which are allowed to move execute their step. 

(5) D is increased by all moving particles. 

The above rules have to be applied to all pedestrians at the same time, i.e. we 
use parallel update. 

Note that we do not use here the so called matrix of preference [6] because 
it encodes a direction of preferred motion of the pedestrians, which they usu- 
ally not have at the beginning of an evacuation process. In the following all 
information about the desired direction of motion is obtained from the floor 
fields. 



5 One can as well solve the arising conflicts using the procedure described in [6]. 
It can be shown [13] that conflicts are important for a correct description of the 
dynamics. 
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In [6] it has been shown that simple, local update rules for the pedestrians 
and the two floor fields are sufficient to yield a richness of complex phenomena 
(see sec. 2.4). Obviously this route is superior concerning the computational 
efficiency and even allows for faster-than-real-time simulations of large crowds 
[6,14,15], e.g. in evacuation processes in public buildings. 



2.4 Collective Phenomena 

Pedestrian dynamics exhibits are variety of fascinating collective effects [2,4]: 

• Jamming: At large densities various kinds of jamming phenomena oc- 
cur, e.g. when many people try to leave a large room at the same time 
[4,6,14,10,15,16]. This clogging effect is typical for a bottleneck situation. 
Other types of jamming occur in the case of counterflow where two groups 
of pedestrians mutually block each other [17,18]. This happens typically at 
high densities and when it is not possible to turn around and move back, 
e.g. when the flow of people is large. 

• Lane formation: In counterflow, i.e. two groups of people moving in op- 
posite directions, a kind of spontaneous symmetry breaking occurs. The 
motion of the pedestrians can self-organise in such a way that (dynamically 
varying) lanes are formed where people move in just one direction [9]. In 
this way, strong interactions with oncoming pedestrians are reduced and a 
higher walking speed is possible. 

• Oscillations: In counterflow at bottlenecks, e.g. doors, one can observe 
oscillatory changes of the direction of motion [9] . Once a pedestrian is able 
to pass the bottleneck it becomes easier for others to follow him in the same 
direction until somebody is able to pass (e.g. through a fluctuation) the 
bottleneck in the opposite direction. 

• Patterns at intersections: At intersections various collective patterns of 
motion can be formed [2,4,19]. A typical example are short-lived round- 
abouts [2,4] which make the motion more efficient. Even if these are con- 
nected with small detours the formation of these patterns can be favourable 
since they allow for a "smoother" motion. 

• Trail formation: Although human and animal trails are formed for rather 
different purposes their structures have some similarities [20,21]. Often hu- 
man trails are formed as a short-cut which makes it attractive to leave a 
paved path. Animal trails usually are related to chemotaxis and mark the 
way to food places. 

• Panics: In panic situations many counter-intuitive phenomena (e.g. "faster- 
is-slower" and "freezing-by- heating" effects [22]) can occur. For a thorough 
discussion we refer to [4,10] and references therein. 
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In [6,14] it has been shown that the new CA model described in the previ- 
ous subsection is - despite its simplicity - able to reproduce these observed 
collective effects. This is essential if one intends to use the model for real 
applications, e.g. the optimisation of evacuation procedures. 



3 Evacuation Simulations 

In the following we describe results of simulations of a typical situation, i.e. 
the evacuation of a large room (e.g. in case of fire). At this we focus on the 
influence of the sensitivity parameters kn and ks on the evacuation times in 
order to identify the different classes of behaviour exhibited by the model. As 
we will see interesting collective phenomena between the pedestrians lead to 
a nontrivial dependence of the evacuation times on k D and ks- In sec. 3.2 we 
investigate possible alternative definitions of the dynamical floor field. Sec. 3.3 
is devoted to a simple optimisation problem, namely evacuation from a room 
with two doors [23]. 

3. 1 The impact of sensitivity parameters 

The value of ks, the coupling to the static field, can be viewed as a measure 
of the knowledge of the pedestrians about the location of the exit. A large ks 
implies a motion to the exit on the shortest possible path. For vanishing ks, on 
the other hand, the people will perform a random walk and just find the door 
by chance. So the case ks <C 1 is relevant for processes in dark or smoke-filled 
rooms where people do not have full knowledge about the location of the exit. 

The parameter ko for the coupling to the dynamic field controls the tendency 
to follow the lead of others. A large value of hp implies a strong herding 
behaviour which has been observed in the case of panics [10]. 

We consider a grid of size 63 x 63 sites with an exit of one cell in the middle of 
one wall. The particles are initially distributed randomly and try to leave the 
room. The only information they get is through the floor fields. Fig. 4 shows 
typical stages of the dynamics for an initial particle density of p = 0.3, which 
means 1116 particles. In the middle picture of fig. 4 a half-circle jamming 
configuration in front of the door is easy to spot. A typical feature of the 
dynamics is a radial motion of 'holes' created by particles escaping through 
the door. 

The evacuation times T and their variances a are averaged over 500 samples 
and strongly depend on the sensitivity parameters k D and ks- Fig. 5 (a) shows 
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Fig. 4. Typical stages of the dynamics: (a) beginning evacuation (i = 0); (b) middle 
stages; (c) end stage of evacuation with only a few particles left. 

the evacuation times for fixed sensitivity parameter ko of the dynamic field 
and variable sensitivity parameter k$ of the static field and fig. 5 (b) shows the 
evacuation times for fixed and variable ks- Fig. 6 shows all corresponding 
variances o and relative variances a/T of the evacuation times of fig. 5. The 
averaged evacuation times are always measured in update time steps. With 
the generic time-scale of the model [6,11,12] of about 0.3 sec/timestep it easy 
to translate that into a real time value. 
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Fig. 5. Averaged evacuation times for a large room (fig. 4) with an initial particle 
density of p = 0.3 and decay constant 5 = 0.3: (a) a = 0.1 and fixed k^; (b) 
a = 0.3 and fixed ks- 



Let us first consider the case of kr> = 0, i.e. no coupling to the dynamic field. 
In fig. 5(a) one can see the influence of ks- For ks — ► the pedestrians do 
not sense the strength of the field. Therefore they do not have any guidance 
through the surroundings and perform a pure random walk which leads to a 
maximal value of the evacuation time for ks — 0. For ks — ► oo they have full 
information about the shortest distance to the door and the evacuation time 
converges towards a minimal value. The movement of the particles becomes 
almost deterministic. Therefore ks can be interpreted as some kind of inverse 
temperature for the degree of information about the inanimate surrounding 
of the pedestrians. 
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Fig. 6. Averaged variances a and relative variances a/T of the evacuation times of 
Fig. 5: (a) 5 = 0.3, a = 0.1 and fixed kjj; (b) 5 = 0.3, a = 0.3 and fixed ks- 

In the same way the sensitivity parameter ko of the dynamic field works as an 
inverse temperature for the information about the virtual trace. If ks is turned 
on from zero to infinity, a non-zero value of kjj only means additional noise to 
the pedestrians and evacuation times increase for higher coupling strength to 
ko (see fig. 5(a)). 

Much more interesting is the behaviour for fixed ks (fig. 5(b)). The evacua- 
tion times saturate at maximal values for growing sensitivity parameter kp 
of the dynamic field. The most interesting point is non-monotonic behaviour 
of T{k,£>) with the occurrence of minimal evacuation times for non-vanishing 
small values of the sensitivity parameter kp of the dynamic field. Therefore a 
small interaction with the dynamic field, which is proportional to the velocity 
density of the particles, represents some sort of minimal intelligence of the 
pedestrians. They are able to detect regions of higher local flow and minimise 
their waiting times. This effect is most pronounced for intermediate coupling 
(ks ?z 0.4) to the static field S, but very weak for strong coupling to S (e.g. 
ks ~ 10), where the evacuation is dominated by the orientation to the inan- 
imate surrounding. It vanishes again for very weak coupling to S (ks = 0.1, 
not shown in fig. 5(b)) where the movement of the pedestrians is similar to a 
random walk. 

If the coupling to the dynamic field is further increased, the evacuation times 
increase again and saturate at maximal values. The interaction with other 
pedestrians becomes more and more unfavourable similar to the arising of a 
panic situation. The weaker the coupling to the static field S is, the higher are 
the evacuation times: the particles than have less information of the inanimate 
surrounding (for example they cannot find the way to the door because of 
smoke in a fire situation). 

The variances a of the evacuation times (fig. 6) show qualitatively the same 
behaviour as the averaged time values. For strong coupling to S (k s ~ 10, 
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kn < 2) a tends to zero, which means that the process is nearly deterministic. 
For weak coupling to S the variance a increases strongly and also the relative 
variance ^ becomes rather large for ks — > 0. From a practical point of view, i.e. 
evacuations of real buildings, one has to ask whether it makes sense to specify 
safety just by the average evacuation time alone when (relative) variances can 
become large. For fixed ks and varying coupling ko to the dynamic field the 
variance a also behaves qualitatively as T, i.e. it increases with increasing ku 
(see fig. 6(b)). 

In the following we discuss the influence of the decay and diffusion parameters 
5 and a. In fig. 7 (a) one can see that for k s = 0.4 and a high density of p = 0.3 




Fig. 7. Averaged evacuation times for a large room for an initial particle density of 
p = 0.3: (a) k s = 0.4 and S = 0.3; (b) k s = 0.4 and a = 0.0. 

the effect of minimal evacuation times for non- vanishing small values of ko 
becomes most pronounced in the limit a — > 0. Besides that all evacuation 
times are increased with increasing a. That implies that diffusion of the field 
only means additional noise and no advantageous information for the particles 
due to the fast broadening and dilution of the trace for large a. 

For decay parameter 5 — > the dynamic field persists for a long time and the 
regime of minimal evacuation times is shifted towards smaller hp values. Evac- 
uation times for fixed kr> increase monotonically with 5 (see fig. 7(b)). This 
is a collective effect where the herding behaviour helps to overcome the insuf- 
ficient knowledge about the location of the exit. In contrast, the evacuation 
times in the panic regime (i.e. high ko values) are increased with smaller 5. 
With the field values of D the memory of previous steps grows which increases 
the noise in the system for the particles. Therefore for large values of ^ a 
strong herding behaviour would be unfavourable since it tends to 'confuse' the 
pedestrians which already have a good knowledge about the geometry. This 
reversed behaviour in comparison with the case k^ — > leads to a crossing of 
the curves in fig. 7(b). It is interesting to note that this crossing happens in a 
rather narrow interval around kp ~ 1. 
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The influence of the diffusion parameter a of the dynamic field in the regime 
of minimal evacuation times strongly depends on the particle density p. As in 




Fig. 8. Averaged evacuation times for a large room (Fig. 4) for two densities, variable 
a and k s = k D = 0.4: (a) p = 0.3, and (b) p = 0.003. 

fig. 7(a) in fig. 8(a) one can see that in the high density regime an increasing 
diffusion parameter a always increases the evacuation times for arbitrary decay 
parameter S, since the diffusion of the field increases the noise for the particles. 
In the low density regime one finds the opposite behaviour. Here the diffusion 
of the dynamic field leads to a higher degree of information for the particles 
and to favourably long-ranged interactions between them (fig. 8(b)). For large 
a the evacuation time increases with decreasing 5 whereas for small a one 
finds just the opposite behaviour. This leads to a crossing of the curves in 
fig. 8(a). The main reason is that for larger densities the evacuation time is 
mostly determined by the clogging which occurs at the exit whereas for small 
densities clogging is negligible. Furthermore for small a the traces are rather 
sharp and so a larger 5 just reduces the field strength. In contrast, for large a 
the traces broaden quickly and are thus diluted such that they form an almost 
constant background. 

Thus three main regimes for the behaviour of the particles can be distin- 
guished. For strong coupling to ks and very small coupling to ko we find 
an ordered regime where particles only react to the static floor field and the 
behaviour than is in some sense deterministic. The disordered regime, charac- 
terised by strong coupling to ko and weak coupling to ks, leads to a maximal 
value of the evacuation time. The behaviour here is typical for panic situations, 
e.g. the herding tendency dominates. Between these two regimes an optimal 
regime exists where the combination of interaction with the static and the 
dynamic floor fields minimises the evacuation time. Here the individuals have 
some knowledge about the location of the exits, but through some sort of 
cooperation the evacuation time T can be optimised. 

Fig. 9(a) shows the influence of an increasing particle density for the evacua- 
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tion times in the three main regimes. For very small densities one finds a very 
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Fig. 9. Average evacuation times for a large room: (a) density regime p € [0,0.9]; 
(b) more detailed look at the low density regime p £ [0,0.1]. 

strong increase of the evacuation times with increasing density (fig. 9(b)). If 
there are only very few particles in the system, which are distributed randomly 
over the lattice at t — 0, they reach the door nearly independently from each 
other during the evacuation. Than the evacuation times are nearly propor- 
tional to the largest initial distance to the exit. In this regime diffusion of the 
dynamic field D provides advantageous information to the particles (see fig. 8 
(b)). In the density regime p ~ 0.005 a queue in front of the door begins to 
form that strongly controls the evacuation time which then increases nearly 
linearly with growing density due to clogging. 



3.2 Alternative definitions of the dynamic floor field 



In the following we want to give a brief discussion of the consequences of sev- 
eral variations and extensions of the definition of the dynamic field D and the 
corresponding coupling parameter kz>. First we concentrate on the effects of 
a negative sensitivity parameter hp < 0, corresponding to repulsive interac- 
tions between the individuals, on the evacuation times. Fig. 10(a) shows the 
averaged evacuation times T for a high density of p = 0.3 and weak coupling 
ks = 0.4 to the static field for several diffusion parameter values a (5 = 0.3 
fixed). For all diffusion parameters a the averaged evacuation time is increas- 
ing for decreasing sensitivity parameter kjj < 0. This is obvious since the 
dynamic field D marks the regions of highest flow in the system (e.g. nearby 
the exit), because of its relation to the velocity density, and therefore an avoid- 
ance of this regions should lead to an increased evacuation time. An increasing 
a value weakens the sharpness of D in the high flow regions and so the effect of 
the avoidance of this regions is moderated, i.e. the evacuation times decrease 
for increasing a and kr> < 0. 
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Fig. 10. Averaged evacuation times for variations of D and hp for a density of 
p = 0.3 and a lattice of 63 x 63 sites: (a) kjj G [—10,10] for D proportional to 
the velocity density; (b) hp G [—10, 10] for a field D proportional to the particle 
density and D proportional to the velocity density. 

More subtle consequences arise for the dynamics of the model if the definition 
of the dynamic field D is changed. Up to now we have considered a dynamic 
field which is related to the velocity density in the sense that it is only altered 
by moving particles (see sec. 2.2). As an alternative we investigate a field 
D which is altered by all particles: in each time step each particle increases 
the dynamic field value of its site by one, i.e. D is proportional to the particle 
density. Fig. 10(b) shows the evacuation times for fixed coupling strength kg = 
0.4 to the static field, variable kjj G [—10, 10] and the two diffusion parameters 
a = 0.0 and a = 0.1 for this density field D (S = 0.3 fixed). As a comparison 
fig. 10(b) shows again the corresponding times for the velocity density field 
D (see fig. 10(a)) with a = 0.0. For positive coupling parameter ko > 0, i.e. 
attractive interaction, one finds monotonously increasing evacuation times for 
all diffusion parameters a. The time values for different a begin to branch for 
kr> > 1 such that for kr> > 5 the evacuation time for a diffusion parameter 
a = 0.0 is three times higher than for a = 0.1. The reason for that becomes 
obvious from fig. 11(c) where a two-dimensional visualisation of the density- 
dependent dynamic field D is shown. The darkest shaded sites correspond to 
regions of highest D values. Since D is related to the particle density, regions of 
jamming due to high particle concentration are signified by high D values (e.g. 
a half-circle jamming in front of a door, see fig. 11(a)). This effect is strongest 
for a = 0.0, because a non- vanishing diffusion parameter only smoothes the 
field gradients. Therefore an attractive interaction to such a particle density- 
dependent field leads to very high evacuation times due to clustering and 
herding effects known from panic situations, especially for vanishing diffusion 
off). 

For negative sensitivity parameter k^ < 0, i.e. repulsive interaction, one finds 
a different behaviour. In the region of G ] — 1,0[ the evacuation time 
for the density-dependent field case is minimized. For smaller k^ values the 
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Fig. 11. Snapshot of a simulation and the corresponding dynamic floor fields D (for 
a = 0, 5 = 0.3): (a) typical middle stage of the dynamics of the original model 
and (b) the corresponding dynamical field D related to the velocity density; (c) 
dynamical field D for a modified model with D related to the particle density. 

evacuation times increase again monotonically. This is similar to the effect 
of minimized evacuation times for the velocity density related field D and 
k D e]0, 1[ (see fig. 5(b) and inset of fig. 10(b)). A slight repulsion from regions 
of high density combined with a weak guidance through the static ground 
field S (ks = 0.4) should lead to a larger flow and to smaller evacuation 
times, as well as a slight attraction to regions of high flow (inset of fig. 10(b)). 
A strong repulsive behaviour supersedes the directed walk through S and 
therefore is counterproductive. It leads again to increasing evacuation times. 
For ko < the evacuation times for the particle density dependent field 
D are always smaller than for the velocity density dependent field D, since 
a repulsion from regions of higher particle density should always be more 
favourable than repulsion from regions of higher flow. However, fig. 10 indicates 
that a density dependent field with repulsive interactions behaves in some 
aspects as a velocity density dependent field with attractive interactions. 

3.3 Room with two doors 

As a simple example for safety estimations in architectural planing we inves- 
tigate how evacuation times change if a gap between two doors is increased 
from zero to a maximal value. We consider again an ordinary room with no 
internal structure of grid size 102 x 102. We start with one door of width two 
cells in the middle of one wall. Then the door is split into two doors of one 
cell each separated by gaps ranging from 2 to 98 (fig. 12(b) and (c)). Fig. 12 
shows typical stages of the dynamics for different gap sizes. Averaged evac- 
uation times are measured for all three main regimes introduced in sec. 3.1. 
Depending on the regime, fig. 13(a) shows a strong influence of the size of the 
gap for small and for large gaps. In the ordered regime, the evacuation time 
is almost uninfluenced by the gap size. The reason is the almost deterministic 
motion of the particles which do not interact with each other through the 
dynamic floor field D. In the other two regimes, however, the value of the gap 
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Fig. 12. Typical clogging configurations for a room with two doors and (a) no gap 
between doors, (b) medium gap size, and (c) doors at the boundaries. 
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Fig. 13. (a) Evacuation times for all three regimes as function of the gap between 
doors, (b) Averaged flows for selected gap sizes for the three regimes. 

can have a large influence due to the strong interaction effects through D. For 
both the optimal and the disordered regime one finds minimized evacuation 
times for a wide area of the gap length ranging from about 20 to 80 cells. For 
small gaps the evacuation times increase due to the interactions between the 
pedestrians. For large gaps the presence of the side walls has a negative effect 
(fig. 12(c)) and so again the evacuation times increase. For intermediate gaps 
the crowd of pedestrians will be subdivided into smaller groups (fig. 12(b)), 
leading to more favorable interactions between them. This picture is confirmed 
by the averaged flow values J in fig. 13(b). A gap size of 40 cells between the 
doors leads to maximised flow values in the optimal and disordered regime. 



4 Conclusions 



We have studied simple evacuation processes using a recently introduced 
stochastic cellular automaton for pedestrian dynamics which implements in- 
teractions between individuals using an idea similar to chemotaxis. Due to its 
simplicity the model allows very high simulation speeds and is very well suited 
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for the optimization of evacuation procedures even in complex situations. 

We have focussed on studying evacuation times in a very simple evacuation 
scenario. The main purpose was to elucidate the influence of the various model 
parameters, especially of the coupling strengths k$ and kp to the static and 
dynamic floor fields, to facilitate their interpretation. 

Most important are the coupling parameters ko and ks to the dynamic and 
static floor fields, respectively, ks is a measure for the knowledge of the individ- 
uals about the geometry, especially the location of the exits. For dominating 
coupling to the static field the pedestrians will choose the nearest exit without 
much detour. For ks — > 0, on the other hand, they have no knowledge at all 
and will find the doors just by chance. 

The dynamic floor field D is a measure for the velocity density of the pedestri- 
ans. The parameter hp controls the herding behaviour. For dominating cou- 
pling to the dynamic field the pedestrians have a strong tendency to follow in 
the footsteps of others, e.g. because they hope that others have more knowl- 
edge about the location of exits. Such a behaviour is relevant for panic situa- 
tions where this herding tendency becomes important and has been observed 
empirically [4,10]. In fact, the ratio kc/ks has very similar effects as the panic 
parameter introduced in [10]. 

The investigation of alternative definitions of the dynamic floor field has shown 
that the most realistic results are obtained for the velocity density dependent 
field with repulsive interactions as used in the original definition of the model. 
However, a field related to the particle density with repulsive interactions 
behaves in some aspects very similarly. 

An important result of our investigations is that for achieving optimal evacu- 
ation times a proper combination of herding behaviour and use of knowledge 
about the surrounding is necessary. Then through cooperative behaviour of 
the individuals evacuation times can be minimized. 

We want to emphasize that here only a very simple scenario has been studied 
in order to identify the different regimes of the model. For realistic applications 
a procedure to determine the coupling parameters is needed. This work is in 
progress and results will be published elsewhere [13,7]. In the Appendix one 
possible idea is discussed (see also [6]). Furthermore the effects of disorder (e.g. 
different individuals j having different couplings k^ and k^) are important. 
Another interesting open question concerns the details of the dynamics of 
evacuation processes in panic situations. We have assumed, as is also done 
for the panic parameter in [10], that the coupling constants do not change 
during the process. However, in many realistic situations it might be more 
appropriate to use a panic parameter which increases with time. We leave this 
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important problem for future study. 

It is suprising that the properties of the model investigated here are in many 
respects very similar to the social-force model [9] although the interactions 
are very different. In our approach pedestrians interact with the velocity- 
density through an attractive coupling whereas the social-force model uses a 
repulsive density-density interaction. It would be interesting to get a deeper 
understanding of these similarities. 
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A Construction of the static floor field S 

The values SV,- for the static floor field for sites for the geometries used 
in this paper can be calculated in the following way: The lattice shall be 
surrounded by obstacle cells, which are not traversable, except for a number 
of door cells at {(i Tl ,j Tl ), • • • , (iT k ,jT k )}- The pedestrians can only leave the 
room through these door cells. The explicit values of S in the examples studied 
here are than calculated with a distance metric: 

Sij = min \ max j yj (i Ta - it) 2 + {j Ta - jt) 2 I - y] (i Ta - i) 2 + (j Ta - j) 2 \ ■ 

(A.I) 

This means that the strength of the static floor field depends on the shortest 
distance to an exit, max^^) 

y^T s -ii) 2 + (j Ta -Ji) 2 }, where (i hJl ) runs 
over all cells of the lattice, is the largest distance of any cell to the exit at 
(iT s ,jr s )- This is just a normalisation so that the field values increase with 
decreasing distance \J (Jr s — i) 2 + {jr a — j) 2 to an exit and is zero for the cell 
farthest away from the door. 

This is only one possible representation of the explicit calculation of the static 
floor field S. Alternative constructions of S do not change the qualitative 
results of the investigations, provided that the strength of the static floor field 
is increased in the direction to the exit [7]. The use of a Manhattan metric 
is one possible example for a successful construction of S for more complex 
geometries [13,7]. 
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Once the floor field is specified from the geometry of the problem the cou- 
pling constant k$ has to be determined such that a quantitative agreement 
with experimental data is obtained. Currently we are investigating different 
strategies [13,7]. A simple method would be to determine ks such that the ob- 
served velocity and velocity fluctuations of the motion of a single pedestrian 
are reproduced (see [6]). 
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